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Abstract 

We develop technique concerning statistical analyses of the large-scale structure of 
the universe in the presence of biasing in the structure formation. We formulate the 
diagrammatic method to calculate the correlation functions of the nonlocally biased field 
from a generally non-Gaussian density field. Our method is based on generalized Wiener- 
Hermite expansion of the density field. The present formalism has not only general 
applicability, but also practical significance, too. To show the effectiveness of this method, 
we revisit the problems on biasing that have been considered previously, i.e., various 
approximations for peak statistics and hierarchical structure of correlation functions of 
locally biased field. Further analyses which have not been possible so far on these problems 
can be performed by our formalism. The gravitational evolution of primordial fiuctuation 
or other structure-forming processes can be described by nonlocal biasing, so can be 
treated by our new formalism in principle. 

Subject headings: cosmology: theory — galaxies: clustering — galaxies: formation — 
gravitation — large-scale structure of Universe 
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1 INTRODUCTION 



The statistics of the large-scale structure of the universe is one of the widest window 
through which we can search for the origin and evolution of our universe. Various mod- 
els on generating the primordial fluctuation of the universe predict different statistical 
properties of the universe. For instance, simpler inflationary models naturally predict the 
Gaussian initial fluctuation with scale-invariant primordial power spectrum (Guth & Pi 
f982; Starobinskii f982; Hawking f982; Bardeen, Steinhardt & Turner f983). On the 
other hand, topological defects (see, e.g., Kolb & Turner f990; Vilenkin f985) such as 
global monopoles, cosmic strings, domain walls or textures produce various non- Gaussian 
fluctuation as well as various non-simpler inflationary models (Allen, Grinstein & Wise 
f987; Kofman & Linde f985; Ortolan, Lucchin & Matarrese f989; Salopek & Bond f990; 
Salopek, Bond & Efstathiou f989; Hodges et al. f990; Salopek f992), the explosion model 
(Ikeuchi f98f; Ostriker & Cowie f98f), and so on. Density fluctuation generated by 
various models determines the statistical properties of mass distribution of the present 
universe which are subjected by the subsequent gravitational evolution. 

The mass distribution, however, can not be observed directly. We can directly ob- 
serve the luminous cosmic objects, such as galaxies, clusters of galaxies, QSOs, X-ray 
clusters, etc. Generally, the number density of luminous cosmic objects of some kind is 
not proportional to the mass-density field and is biased tracers of mass. The relation 
between the number density field pL{f ) of the luminous objects and the density field 
Pb{x) of underlying mass has the essential effects on determining the mass distribution 
of the present universe. These two fields can be considered as nonlocal functionals of a 
primordial fiuctuation field. If we can express the relation of those two fields, the relation 
generally should be a functional as Phif) = Fr[Ph{x)\- We assume this type of relation 
exists in this paper and explore the statistical consequences of this relation of biasing. 
This functional depends on the physical process to form the luminous objects and the 
primordial fiuctuation of the universe. 

As an example, in usual biased galaxy formation scenarios (Kaiser 1984; Davis et 
al. 1985; Bardeen et al. 1986; Coles 1989; Lumsden, Heavens & Peacock 1989), the galax- 
ies form at high density peaks of primordial fiuctuation. Adopting this hypothesis, the 
standard cold dark matter (CDM) model can explain the observations of two-point corre- 
lation function (e.g., Davis & Peebles 1983) on small scales. However, recent observations 
as Automatic Plate Measuring Facility (APM) survey (Maddox et al. 1990), Queen Mary 
and West field College Durham, Oxford and Toronto (QDOT) redshift survey (Ffstathiou 
et al. 1990; Saunders, Rowan- Robinson & Lawrence 1992; Moore et al. 1992) suggest that 
the clustering of galaxies on large-scales ^ fO/i~^Mpc is stronger than that of standard 
CDM model. It is pointed out that this discrepancy can be reconciled by introducing 
nonlocally biasing mechanisms as quasor-modulated galaxy formation (Babul & White 
1991) or cooperative galaxy formation (Bower et al. 1993) rather than the simple biasing 
by peaks. 

Though the biasing by peaks is a simple model to investigate, more complicated bi- 
asing is certainly present in reality. The real bias is determined by gravitational and 
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hydrodynamical mechanisms (Cen & Ostriker 1992, 1993), but the right form of nonlocal 
bias is not known yet. 

The local bias, represented by a function pL{f ) = f{pb{f )) instead of a functional, is 
widely explored so far. Kaiser (1984) showed that the enhancement of two-point corre- 
lation function of clusters relative to that of galaxies (see, e.g., Bahcall & Soneira 1983; 
Klypin & Kopylov 1983; Hauser & Peebles 1973; Postman, Huchra & Geller 1992) can 
be explained by local biasing between galaxy distribution and cluster distribution. He 
adopted a model that the biased objects are homogeneously formed in the region where 
the value of underlying field 6}j{r) is larger than the threshold v in units of rms a of 
underlying field, and then calculated the two-point correlation function of the biased ob- 
jects in high and large- separation limit for random Gaussian underlying field. In this 
model, the biasing is local and represented by a step function as pB{f ) = d{8i{r) — per). 
The generalization of this result to A^-point correlation functions or to the expressions 
with loosened restriction on p and/or separation (Politzer & Wise 1984; Jensen & Szalay 
1986) was performed. Further generalizations to the case in which the underlying field 
has non-Gaussianity (Fry 1986; Grinstein & Wise 1986; Matarrese, Lucchin & Bonometto 
1986), and/or the case of general form of local biasing function which is not necessarily the 
step function (Szalay 1988; Borgani & Bonometto 1989, 1990). In these generalizations, 
however, the biasing is restricted to be local. 

The biasing by peaks does not fall under the category of local bias. The average peak 
number density for Gaussian random field was analytically calculated (Doroshkevich 1970; 
Adler 1981; Bardeen et al. 1986). Some approximations for correlation functions of peaks 
are known (Bardeen et al. 1986; Otto, Politzer & Wise 1986; Cline et al. 1987). Average 
peak number density for generally non- Gaussian random fields for high-threshold limit is 
also known (Catelan, Lucchin & Matarrese 1988). 

In this paper, we generalize those previous works from an unified point of view. The 
mathematical tools for investigating the relation between biasing and statistics in gen- 
eral situations are presented. As the general framework, the mean density and A^-point 
correlation functions for biased fields are expressed as series expansions. The underlying 
field is not restricted to be Gaussian, and the biasing is not restricted to be local. The 
diagrammatic methods to give each terms in the series expansion are introduced. In §2, 
the generalized Wiener- Hermite functionals which play essential roles in our analysis are 
introduced. Average value of any functional of a field is expanded by Wiener-Hermite 
functionals. In §3, the diagrammatic method to calculate the mean density and A^-point 
correlation functions of biased field is derived. More compact expressions or methods of 
calculation is given when the biasing is local or semi-local, which is defined below. The 
Fourier-space version of our methods are also explained. Our methods are applied to 
several theoretical problems in §4. The first application is not relevant to bias, but the 
relation to Fdgeworth expansion is discussed. The Fdgeworth expansion is recently used 
in analyses of astrophysical density fields. The second application is on peak theory. We 
show that our methods improve the techniques which is available so far for peak statistics. 
The next application is on hierarchical relations of correlation functions in the presence 
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of local biasing which is investigated recently by Fry & Gaztanaga (1993; and references 
therein). We show our method complements their results. The last application is on the 
gravitational nonlinear evolution as a nonlocal biasing. We apply our methods to the 
calculation of three-point correlation function induced by gravity. Finally, we present the 
conclusions in §5. Appendices are devoted to the details of arguments and calculations. 

The main purpose of this paper is to introduce a general framework to calculate the 
statistics of nonlocal biasing from non-Gaussian fields. In this paper, applications are 
restricted to the previously analyzed problems. Our methods are fairly general, and it is 
applied not only to the problem of biasing but also to the statistical analyses of large-scale 
structure of the universe. 



2 NON-GAUSSIAN RANDOM FLUCTUATIONS 



2.1 Correlation statistics 

We introduce first a homogeneous random field a[x) with zero mean in three-dimensional 
space. For astrophysical applications, this random field corresponds to, e.g., the smoothed 
density fiuctuations, 

Snix) = J d'yWni\x-y\)Siy), (2.1) 

or normalized density fiuctuation, 6r[x) / {S^y^'^ ^ etc., depending on specific applications. 
In the applications in this paper, a is always regarded as normalized density fiuctuation. 
In the above notations, 6 = p/ p — 1 is the density contrast and Wr is the window function 
which cut the high frequency component higher than 1/R. Two popular forms are the 
Gaussian window: 

>»t'(..)=^5^exp(-5^1), (2.2, 

and the top-hat window: 

wr\^) = ^AR-^)- (2-3) 

We review some correlation statistics for a general field a in this subsection. 

The probability distribution for a field a is specified by a functional ^[a]. Averaging 
(• • •) is represented by a functional integration as 



(•••) = J[da]...V[a]. (2.4) 

The mean value of A^-product of the field is called A^-th moment: 

pS^\xi,. . . , Xn) = {a{xi) ■ ■ ■ a{xN)). (2.5) 
These moments are generated by the following generating functional: 



Z[J] = j [da]V[a] exp i j d^xJ{x)a{x 



(2.6) 
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through the functional differentiation: 



6J[xi) ■ ■ ■ 6J[x 



N 



(2.7) 



j=o 



We usuaUy use reduced moments rather than mere moments. The A^-th reduced moment 
is defined by 

V'^^^a?!, . . . , Xn) = {a{Xi) . . . a(a;jv))connected, (2.8) 

where (• • •)connected dcuotcs "conuected part" of the moment (Bertschinger f992) removing 
the disconnected parts from the moment. The reduced moments are called semi- invariants 
or cumulants in probability theory and related to the connected Green's function in quan- 
tum field theory. If a is identified with the density contrast (5, the A^-point correlation 
function in the continuum limit in cosmology (e.g., Peebles 1980) is identical to the 
A^-th reduced moment, neglecting shot noise effect. When the statistics of a field a is 
Gaussian, the reduced moments larger than or equal to third order are all vanish. 

According to the celebrated "cumulant expansion theorem" (e.g.. Ma 1985), reduced 
moments are generated by logarithm of moment-generating functional Z[J]: 



\N sN 



6'^ In Z 



Xn 



6J[xi) ■ ■ ■ 6J[x 



N 



(2.9) 



j=o 



where r/j^^^ is set to zero formally. Thus, 



Z[J] = exp 



oo 



E 

.JV=2 



A^! 



1? ■ 



,XN)Jixi) ■ ■ ■ J{X 



N 



(2.101 



provide a very useful way to tackle non-Gaussian random fields because it relates the 
probability distribution functional ^[a] and the reduced moments through the equation 
(2.6). 



2.2 Generalized Wiener-Hermite functionals 

We define generalized Wiener-Hermite functionals which play essential roles in our 
analyses below. For a given random field a(a;), the generalized Wiener-Hermite function- 
als are defined by 

(«!,...,«„) = exp - j d^xd^ya{x)tlj-^x,y)a{y) 



6a{xi) ■ ■ ■ 6a{Xf 



exp 



d^xd^ya{x)'ijj ^{x,y)a{y) 



(2.11] 



where t/j ^ denotes an infinite dimensional inverse matrix of t/j defined formally by 

J d^yr\x,y)i,{y,z) = 6%x - z), (2.12) 
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and t/j is an arbitrary function of two points in space in general situations. In the following, 
we designate t/j in the above definition of generalized Wiener-Hermite functionals as 2nd 
reduced moment r/j^^^ defined by equation (2.8) with N = 2. The generalized Wiener- 
Hermite functionals are reduced to the Wiener-Hermite functionals in the special case 
t/j[x^y) = 6^[x — y). The following functionals defined by infinite dimensional linear 
combinations of equation (2.11) as 

'H(^){xi, ...,«„) = y" J^j/i • • • d^y^ip{xi,y^) ■ ■ ■ ip{x^,y^)n^"'\y^, . . . , ?/„), (2.13) 

are convenient for our purposes below. Note that we distinguish two classes of functionals 
(2.11) and (2.13) by the place of (m). The latter functionals are also called generalized 
Wiener-Hermite functionals. The generalized Wiener-Hermite functionals defined above 
are the natural generalization of the generalized Wiener-Hermite polynomials (e.g., Appel 
& de Feriet 1926) which have finite degrees of freedom rather than three dimensional 
continuum degree of freedom. The first several forms of 7i(m) are 

^(0) = 1, (2.14) 

n(,){x) = a{x), (2.15) 

^(2)(«i, X2) = a{xi)a{x2) - V'(«i, «2), (2.16) 
^(3)(«i,«2,«3) = a{xi)a{x2)a{x3) 

- [ijj{xi,X2)a{x3) + V(«i, «3)a(«2) + V'(«2, X3)a{xi)] , (2.17) 
^(4)(«i,«2,«3,«4) = a{xi)a{x2)a{x3)a{x^) 

- [tlj{xi, X2)a{x3)a{x4) + V(«i, «3)a(«2)a(«4) + V'(«i, «4)a(«2)a(«3) 
+ tp{x2, X3)a{xi)a{x4) + V'(«2, X4)a{xi)a{x3) + V'(«3, X4)a{xi)a{x2)] 

+ [?/'(«!, a;2)V'(«3, X4) + V'(a;i, a;3)V'(a;2, X4) + V'(a;i, a;4)V'(a;2, X3)] . (2.18) 

The generalized Wiener-Hermite functionals Ti^™^ ^(m) are symmetric about its argument 
a?!, . . . , Xfn and are regarded as functionals of a field a[x). Under linear transformation 
of a field a: 

a'{x)= [ d^yM{x,y)a{y), (2.19) 



the generalized Wiener-Hermite functionals transform as follows: 

T~t ^ ("Kl , . . . , Xjyi ) 

= j d^y4---d^y^M-\x4,y4)---M-\x^,y.^)H^^\y4,...,y.^), (2.20) 

= j d^j4 ■ ■ ■ d^j^M{x4,y4) ■ ■ ■ M{x^, y.^)n(^){y4, . . . , j/„), (2.21) 

where is an infinite dimensional inverse matrix of M: 

J d^yM-^x,y)Miy,z) = S^ix - z). (2.22) 
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Thus, we can regard ?-{(™) as infinite dimensional covariant tensors and Ti^m) as infinite 
dimensional contravariant tensors. In physical applications, the linear transformation 
(2.19) is considered as smoothing of a field a if Af is a low-pass filter or is considered as 
the Fourier transformation if M[x^y) = exp[—ix ■ y). The following recursion relation. 



a Xr 



d^ytlj{x^+i,y) 



8 



6a{y) 



n 



(m) 



X 



1? 



Xrf 



n 



X 



m + 1 J 1 



(2.23) 



are helpful properties of the functionals. 



2.3 Expectation values 

The averaging (2.4) for a non-Gaussian random field is reduced to the averaging for a 
Gaussian random field using the generalized Wiener- Hermite functionals as we shall see 
below. First of all, with a given distribution ^[a], we associate a Gaussian distribution 
"PGict] constructed by the 2nd reduced moment of the original distribution ^[a]: 



PqW] = iVaexp 



cPxcPya{x)ip ^{x,y)a{y] 



where Na is a formal normalization constant: 



Nc 



[da] exp 



d^xd^ya{x)il) ^ 



x,y)a[y, 



(2.24) 



(2.25) 



The averaging by this Gaussian distribution "PGict] is denoted as (• • •)g distinguished from 
the averaging by the original non- Gaussian field: (• • •). 

The relation between the non- Gaussian distribution ^[a] and the Gaussian distribu- 
tion Via] is 



Via] = exp 



.JV=3 



d^x^---d^xml^^^\xi 



8 



8 



8a{xi] 



8a{x 



N 



PgW]. 



(2.26) 

This relation is proven as follows: the generating functional Z[J] is (infinite dimensional) 
Fourier transform of distribution functional ^[a]. Thus the inverse Fourier transformation 
gives the expression of ^[a] in terms of reduced moments because of equation (2.10). 
Substituting i8/8a[x) for J(a;), and performing Gaussian integration, we finally arrive at 
the relation (2.26). Fxpanding the exponential of equation (2.26) and using the definition 
of the generalized Wiener-Hermite functionals, we can express non-Gaussian averaging of 
arbitrary functional F[a] by Gaussian averaging as follows: 



{F[a 



: {F[a])o 

CO 1 CO 



E 

'TTi — 



1 m =3 



3 

(1) 



X 



(ni) 



1=1 k=l 
(1) 



, Xl 



1 = 1 

m J \ 



a 



G 



X 



(2.27) 
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This equation is an expansion of a non- Gaussian averaging by higher order reduced mo- 
ments. The nontrivial part in the equation (2.27) is the factor of the form, {TiF)^. In 
Appendix A, the methods to calculate the Gaussian averaging of the product of general- 
ized Wiener-Hermite functionals, (n^)G are developed. Therefore, if a functional F[a] 
is expressed by the sum of the product of generalized Wiener-Hermite functionals, the 
right-hand-side of equation (2.27) can be calculated as performed in the next section. 

3 DIAGRAMMATIC CALCULUS OF TV-POINT CORRELATION 

FUNCTIONS 

3.1 Generalized Wiener-Hermite expansion and diagrammatic methods 

In this paper, we are interested in the correlations of a density field /Ob(^7 [ct]) which 
appears by some biasing mechanism from a field a[x). We formally do not restrict 
the form of biasing in this subsection and consider a given, generally nonlocal, biasing 
functional /Ob(^', [«]) which is both a function of a point r and a functional of a field a[x). 
It is natural to expand this functional by the generalized Wiener-Hermite functionals as 
follows: 



OO 1 . 

n ml J 

m = U 



tV ..... 



(3.1) 



where the kernel A^(™) is given by 

A^(")(r-a;W,...,r-a;(")) = {n^^\x^^\ . . . ,x^^^)pB{r,[a]))^ (3.2) 



(5a(a;(i)) • • •(5a(a;(™)) 



(3.3) 



The equation (3.2) is derived from the orthogonality of the generalized Wiener-Hermite 
functionals (A. 20). We assume that the expansion (3.1) exists for the nonlocal biasing of 
our interest. Using equation (2.27), the A^-th moment of a biased density field is given 
by somewhat complicated expression: 

Pjv(ri, . . . , rjv) = {pB{ri, [a]) ■ ■ ■ /jb(»'jv, [a])) (3.4) 

N m, N 



OO OO .. 1\ mi ^'1 



mi=0 mjv=0 ?'=1 a = l «'=1 



/ \ OO -1 OO OO 

n^™,(a^f\...,a^S"'M + E Ae--- E - 

\i = \ I G m = l n\=i nm=-i ^ 

m Til Tfi 

1=1 k=i 1=1 



/lib "J I III, 
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/ N 



.(3.5) 



This expression can be calculated making use of the methods described in Appendix 
A. The diagrammatic method in calculating Gaussian expectation value of the product 
of generalized Wiener-Hermite functionals is translated to the diagrammatic method in 
calculating moments (3.5) of p^. The result is summarized as the following procedures: 

Draw N open circles corresponding to ri, . . . , rjv which we call external points. 



A-i 
A-h 



A-hi 



A-iv 



A-v 



Consider one of possible graphs with lines between external points or between ex- 
ternal points and vertices. Vertices are the parts where three or more ends of lines 
are gathered. Do not connect two vertices directly with lines. Do not draw any line 
whose two ends are placed at the same external point or the same vertex. 

Associate labels like a;,?/, etc. with ends of lines placed at external points (not at 
vertices). 

Apply the correspondence of Figure 1 for the external points, lines between external 
points and vertices. Make the product of these factors of the graph. 

Integrate with respect to all the labels and multiply the following statistical factor: 

^ (3.6) 



ml Yi «endJ 
ends 

where m is the number of vertices and agnds is the number of lines whose ends are 
identical. This statistical factor is for the case we distinguish the vertices in the 
graph, i.e., we distinguish the graphs like Figure 13 described in Appendix B (see 
Appendix B for more clear understanding of this situation). 

A-vi) Sum up values obtained from all the possible graphs. 

The readers are invited to calculate several examples both by equation (3.5) and by the 
above diagrammatic rules. The Appendix B is devoted to the explanation of the statistical 
factor (3.6). 

This diagrammatic rules enable us to easily calculate lower order terms about ?/'(™) 
or A"("). Although we are not guaranteed the sufficient convergence of this diagram- 
matic expansion in the general situation, there are many cases in astrophysics where this 
expansion is useful as we shall see in §4. 

The reduced A^-point correlation function of the biased density field pB- 

Clf Vl, • • • . ^iv) = i {pBiri, [a]) ■ ■ ■ pBivN, [a])) connected (3-7) 

Pb 
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is more easily evaluated in this diagrammatic expansion because the contribution to the 
reduced moment in equation (3.7) comes from connected graphs in our approach. This 
especially simplify the calculation of higher order correlation functions of the biased field. 
Thus, the above rules are applicable to the calculation of reduced A^-point correlation 
functions substituting the following rule for A-vi). 

A-vi) Sum up values obtained from all the connected graphs and multiply {Pi)~^ . 



3.2 Local and semi-nonlocal biases 

Here we consider special cases of the above general consideration, i.e., the case biasing 
is local or semi-nonlocal. 

The local bias is defined so that the biasing functional /Ob(^7 M) is a mere function of 
the underlying field a at the position r : 

pB = pB{a{r)). (3.8) 

In this case, the kernels [eq. (3.3)] have the forms, 

A^(")(r - a;W, . . . , r - «(")) = S^{r - a;W) • • • S^{r - a;("))i?„, (3.9) 

where 

If the field a is normalized so that (a^) = 1, are coefficients of the Hermite expansion 

of the local biasing function: 

OO 

PBia) = Y] —^H^{a); (3.11) 

m = U 

1 2 

i?„ = — / dae-" l^pB{a)H^{a). (3.12) 
The definition of Hermite polynomials and its orthogonality relations are, 

HUa) = e"'"[-j^^ e-"'l\ (3.13) 

and 

1 /"°° 2 

Jxe"^ '^Hm{x)Hn{x) = ml S^n, (3.14) 



respectively. If the field a is not normalized and (a^) = cr^ 7^ 1, the corresponding 
expansion is 



PB[a) 



J2 ^(T^H^{a/a); (3.15) 



m = U 



Rm = I ~ dae-"'^('^"^pBia)H^ia/a). (3.16) 
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The special form of equation (3.9) reduce the diagrammatic rules A-i) ~ A-vi) above to 
the following: 

B-i) Same as A-i) 
B-ii) Same as A-ii) 

B-iii) Apply the correspondence of Figure 2 and make the product of these factors of the 
graph. 

B-iv) Multiply the statistical factor (3.6). 

B-v) For Pjv, sum up values obtained from all the possible graphs. For reduced A^-point 
correlation function, sum up values obtained from all the connected graphs and 
muhiply (Pi)-^. 

The above rules for Pjv could be expressed by the following equation: 

oo oo oo oo oo o r> 

PMrr,...,r^)=j: ^ •.• ^ •.• ^ •.• ^ ""'^l "J^'''^''' 

m = Oran=0 raijv=0 nml=0 n^jf=0 J~J ]^n,-fc! 

j = l k=l 

u;(^'='^")(l, . . . , 1, . . . , A^, . . . , A^) • • • u;(^'='^'"'=)(l, . . . , 1, . . . , A^, . . . , A^). (3.17) 

nil niN ™ml nmN 



where 



w 



(°) = o 



u;W(z) = 0, 

w^^Hh, . . . , z„) = V^'^H^n ■■■,rj (n > 3). (3.18) 



Identifying the underlying field a with normalized density fiuctuation 6ji[x) / {S^y^'^ ^ 
this formula is an alternative expression to one obtained by Borgani & Bonometto (1989, 
1990) who adopted the different approach to treat the local biasing function pB in eval- 
uating the same quantity of equation (3.17). Our expression seems to be simpler than 
theirs. 

When the statistics of an underlying field is random Gaussian, equation (3.17) reduces 

to 

Piv(n,...,r^.) = f: E n ^'^^'^'^r^' -n^^^.. (3.19) 

m = {mki\m} k<l ^^l' r=l 

In the above equation, 1 < A;, / < A^, 



r-l N 

m^ = ^rnkr+ E ^-rk (3.20) 

k=l k=r + l 
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and Y,{mki\m} means the sum over mki's[k < I) which are non-negative integers with a 
constraint, J2k<i''^ki = We can see, using multinomial expansion theorem, that the 
equation (3.19) is equivalent to the formula obtained by Szalay (1988). 
When the local biasing function pB is given by sharp clipping: 



C9{ 



a 



(3.2i: 



then the coefficients are given by 

r c 

/?(sc) _ 



m 



0) 



2''''' [7=2) 
^ H^_,{v)e-^'l^ (m>l) 



(3.22) 



where 



erfc(a;) = 1 — erf (x) 



TT Jx 



'dt 



(3.23) 



is a complementary error function. In this case, equation (3.17) reduces to the formula 
obtained by Matarrese et al. (1986). If we, furthermore, constrain the statistics of under- 
lying field to be Gaussian, and local biasing function to be the sharp clipping, we can find 
equation (3.19) reduces to the the formula obtained by Jensen and Szalay (1986). 

High-;/ limit of equation (3.17) for sharp-clipping biasing function (3.21) enables us to 
sum up infinite series of graphs and results in 



c 



N 



N 



exp 
G 



00 1 
m = 2 ni,...,njv>0 ^1- ^iV- 



u;(™)(l,...,l,...,iV,...,iV) 



"1 H h" w = "i 



JV 



^27rz/e'' 



where 



Zn (Ji 




(3.24) 
(3.25) 

(3.26) 



is a A^-point moment generating function Z]\j{ Ji) in which v is substituted for all Ji. This 
formula is equivalent to the one obtained by Matarrese et al. (1986). The N = 2 case of 
equation (3.24) has the same form obtained by Grinstein & Wise (1986). 

The local bias is not relevant to the biasing through density peaks (Davis et al. 1985; 
Bardeen et al. 1986; Coles 1989; Lumsden et al. 1989) because density peaks are described 
by not only the mere value of the underlying field but also spatial derivatives of it (see 
§4 for details). Thus, we next consider the case that the biasing functional /JB(^',[ct]) 
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depends on a[r) and finite-order derivatives of it at r, generaiiy. We denote tlie spatiai 
derivatives of a underlying field a as a^: 



/ \ 

da d^a 

Q »Q < j), • • • 

\ 1 ^ ^ ' / 

\ 3 6 / 



(3.27) 



where // = 0, 1, 2, . . .. This notation enable us to write the semi-nonlocal biasing functional 
as 

PB{r,[a]) = pB{a,{r)). (3.28) 
In this case, kernels (3.3) have the forms, 

K(-\r - . . . , r - a.(™)) = E " " " E -^..(^ " ^^'^) ■■■^.Jr- a.(™))i?,,...,„, (3.29) 



/til Aim 



where (5^ denotes the derivatives of Dirac's delta-function in the notation of equation 
(3.27), and 

/ d d \ 
Rf^i-f^m = { ^ Pb{o:^,)) . (3.30) 



The special form of equation (3.29) reduce the diagrammatic rules A-i) ~ A-vi) above to 
the following: 

C-i) Same as A-i) 
C-ii) Same as A-ii) 

C-iii) Associate labels like etc. with ends of lines placed at external points. 

C-iv) Apply the correspondence of Figure 3 and make the product of these factors of the 
graph. 

C-v) Perform sums with respect to all the labels and multiply the statistical factor (3.6). 

C-vi) For Pjv, sum up values obtained from all the possible graphs. For reduced A^-point 
correlation function, sum up values obtained from all the connected graphs and 
muhiply (Pi)-^. 

In Figure 3, V'^^H^" '''j) denotes the derivatives of r/j^^^ and ji indicates the derivatives with 
respect to Vi according to equation (3.27), v indicates the derivatives with respect to 
Similar definition is applied to V'^"1..a(^«5 •••?'''*:)• 
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3.3 Fourier transforms 

The biasing may be expressed in Fourier space rather than a coordinate space. For 
example, the transler lunctions T[k) of various cosmological models is a kind of biasing 
in our context through Fourier space. It may be useful for the analyses in the future to 
give the diagrammatic rules in Fourier space. 

The Fourier transformation is a kind of the following linear transformations of a field: 



a'{k) = j (PxM{k,x)a{x) (3.31) 

PBip.W]) = Jd\Liv,r)pBir,[a]), (3.32) 

where M and L correspond to the coefficients of the linear transformations. Fourier 
transformations are achieved by setting 

M{k, x) = exp{—ik ■ x), L{p,r) = exp{—ip ■ r), (3.33) 

but we use general notations for a while. Because of the property (2.20), (2.21) of the 
generahzed Wiener-Hermite functionals, the form of equation (3.5) is unchanged under 
hnear transformations (3.32), if we demand the following transformation properties: 

V''('^)(fci, . . . , fc„) = I d^x^ ■ ■ ■ d^x^M{k^, x^)--- M{k^, a;„)V'(")(«i, • • • , Xn). (3.34) 
/r(")(p; fc(^), . . . , fc(")) = J d\ J d^x^^^ ■ ■ ■ d^^'^^^Lip, r) 

X M-\x^^\ fcW) • • • M-^(a;("), fc("))A^^"V " • • • , ^ " «^"^)- (3-35) 
The equation (3.35) has the another expression, 

A-(™)(p;fcW,fc(™)) = / '^>b(p,H) \ (3_3g^ 

' ' \Sa'{k^^^)---Sa'{k^'''^)/^ ^ ^ 

because of equation (3.3). These transformed quantities and A^'(™) are equaUy apphed 
to the diagrammatic rules i) ~ vi) instead of the original quantities. This gives a way 
to calculate the correlations of a biased field expressed by linearly transformed quantities 
(3.32). 

Now we restrict the argument to the Fourier transformations (3.33). In this case, from 
the inverse Fourier transformation, M~^[x^k) = exp[ik ■ x)/[2Try. Then transformed 
quantities are given by 

V''('^)(fci, . . . , fc.) = V^^-^H^i, . . . , fc.) = {27rfS'{k^ + ■■■ + k^)i'^\k^, . . . , fc„_i), (3.37) 

and 

A^'(")(p; fcW, . . . , fc(")) = (27r)3,53(fcW + . • .+fc(™)-p)^_ir(™)(fcW, . . . , fc(™)), (3.38) 

(27r )-^™ 
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where and A^(™) are the Fourier transformations of and K^"^\ respectively. Be- 
cause of the homogeneity of a field a, the last expression of equation (3.37) is possible 
and gives a definition of the multi-spectrum p^") of a field a. In the case n = 2, p^^^ is 
called power-spectrum of the field. The Fourier transform of A^(™) is given by 

. . . , . ,-2.,- §-Mv, ISl)^^ . (3.39, 

from the equations (3.36) and (3.38). This form is convenient if the biasing functional is 
expressed in Fourier space, noting that the Gaussian averaging in equation (3.39) is given 
by 

/ |5(fc)p\ 2\a{k)\d\a{k)\ dO^ 
''''^ V Vp{k) ) Vp{k) ~2^_ 



{■■■h = i--- n 

fceuhs 



(3.40) 



in Fourier space, where uhs means a upper half fc-space, 6/^ is a phase of a[k) so that 
a[k) = |5(fc) I exp(z^^). The power spectrum p(A;) is defined by y~^(|5(fc)p), where V 
is the total volume of the space. Substituting equations (3.37) and (3.38) to our rules, 
moments Pjv or correlations i^g^^ in Fourier space are obtained. The multi-spectrum 
Pb^\pi7 ■ ■ ■ ^Pn-i) which is defined by 

ii^\pi, ...,Pn) = {2^f^\Pi + ■■■+ Pn)p\^'\pi,- ■ ■ ,Pn-i), (3.41) 

is more popular than i^g^^ 

The procedures to obtain Pjv are summarized as the following rules: 
D-i) Draw N open circles corresponding to p^^. . . ^pjy which we call external points. 

D-ii) Same as ii) 

D-iii) Define orientations to all lines and associate labels like fc,Z, etc. to all lines. The 
orientations of a line k is regarded as a direction of a current k. The sum of currents 
fiowing into a vertex should vanish and the sum of currents fiowing from a external 
point i should be p- (these rules are due to the (5-function in equations (3.37) and 
(3.38)). 

D-iv) Apply the correspondence of Figure 4 for the external points, lines between external 
points and vertices. Make the product of these factors of the graph. 

D-v) Perform integration / d^^kj {2ttY for each label k and multiply the statistical factor 
(3.6). 

D-vi) Sum up values obtained from all the possible graphs and multiply {2ttY8'^{p-^ + • • • + 
Pn)- 
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From the value of Pi which can be evaluated by the above rules, the mean density of 
biased field Pi = pB is given by 

Pi{p) = {27rf6%p)pB. (3.42) 

Using the rules for obtaining the multi-spectrum pg^^ of biased fields are as follows: 

E-i) Draw N open circles corresponding to Pi^. . . ,/>jv-i Pn — ~{Pi + ' ' ' + Pn-i) 
which we call external points. 

E-ii) Same as ii) 

E-iii) Define orientations to all lines and associate labels like fc,Z, etc. to all lines. The 
sum of currents fiowing in a vertex should vanish and the sum of currents fiowing 
from a external point i should he p-. 

E-iv) Apply the correspondence of Eigure 4 for the external points, lines between external 
points and vertices. Make the product of these factors of the graph. 

E-v) Same as D-v) 

E-vi) Sum up values obtained from all the connected graphs and multiply (/jb)~^- 



4 THEORETICAL APPLICATIONS 



4.1 The Edgeworth expansion of density probability distribution function 

Recently, the Edgeworth expansion of probability distribution function P{S) of cosmo- 
logical density contrast 6 is suggested to be useful in cosmology (Scherrer & Bertschinger 
1991; Juszkiewicz et al. 1994; Bernardeau & Kofman 1994). This expansion is suitable to 
obtain approximately density probability distribution function when restricted number of 
cumulants are known. The expansion is (Bernardeau & Kofman 1994) 



V27r(72 



S3S4 



120 ^ ' 144 



1296 



+ 



(4.1: 



where a = \J (S^) is rms of density contrast, ly = S / a is a normalized density contrast, and 
Sn = ('^")c/('^^)"~^ are normalized cumulants. We can see the coefficients of Edgeworth 
expansion at arbitrary order are very easily calculated in our diagrammatic methods as 
follows. 
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Density probability distribution lunction expressed as 

P{6) = {6n{S - (4.2) 
is appropriate to our purpose. Setting pg = 6£)[acr — v(t), equation (3.12) reduces to 

i?™ = -^L=e-''"I^HUv). (4.3) 

Because tp^^^r ^ ■ ■ ■ ^f) = Sn(y^~'^ , a vertex with n-lines corresponds to order cr""^. Thus, 
diagrams of Figure 5 (a), (b), (c), (d), . . . give all the coefficients of order (7°, cr^, cr^, cr"^, 
. . ., respectively. Applying rules B-i) ~ B-iv), the expansion (4.1) is easily obtained. This 
example shows that our diagrammatic method can save the labor of calculation. 



4.2 Density peaks of Gaussian random fields 

We focus on the statistics of density peaks in this subsection. Bardeen et al. (1986) 
investigated the statistics of density peaks of random Gaussian fields identifying the den- 
sity peaks of primordial fiuctuation with the sites for galaxy or cluster formation. They 
gave practical approximations to the A^-point correlation functions of the density peaks. 
We show, in the following, that our method applied to density peaks of random Gaussian 
field provides the improved approximations to the correlation functions along the strategy 
of theirs. 

In this subsection, the field a is identified with normalized density fiuctuation 6ji[x) / . 
The number density of peaks of the field a greater than a threshold p is given by (Bardeen 
et al. 1986) 

pAr) = OiMr) - ^)'^' {dMr)) {-i f det {cWMr)) (Ai(r)) 6 {\,{r)) 9 {Xs{r)) , (4.4) 

which is a semi-nonlocal bias in our context. In the above, Ai, A2, A3 are eigenvalues of 
a matrix —(didja). The approximation, which is common to Bardeen et al. (1986), is to 
neglect all the derivatives of the two-point correlation r/j^^^. Other correlations ?/'(^) [N > 
3) do not exist because of the assumption of Gaussianity of underlying fiuctuation. As 
an illustration, let us consider the case the power spectrum of background fiuctuation has 
the power-law form with spectral index n. In this case, the two-point correlation tjj^'^^r) 
falls off as |r I "^^ and the m-th derivative of the function falls off as |7»|-("+3+™). When 
n < —2, all the derivatives of two-point correlation function can be neglected even if 
~ 1. 

Adopting the above approximation, the rules C-i) ~ C-vi) for the semi-nonlocal bias 
of random Gaussian underlying field are reduced to the rules A-i) ~ A-vi) for the local 
bias of random Gaussian field with 
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Note that this simphcity does not hold in non- Gaussian underlying field. Further assump- 
tions for higher-order correlations which can hardly be justified are needed to have this 
simplicity for general non- Gaussian field (see the next section). 

For the density peaks, equation (4.5) can be calculated using equation (A. 18) of 
Bardeen et al. (1986) and results in 



where 



da , 



9(a-iy) 



\ dv 



m — 1 



(m = 0), 



1 





{x — x^f 


exp 


[ 2(1-7^)J 



"^"^-^^"^^ 



+ 



erf 



2 

57r 



-|- erf 



2/ 2 



'31x2 8 
5 



and erf (x) is an error function: 



f(. 



erix 



TT Jo 



'dt. 



Spectral parameters cTj, 7, i?* are defined by 



27r 



(7^ 



O'20'O 



where P(A;) is a power spectrum: 

m 



(7 



2 / J3 



jfc-a; / (2) 



(4.6) 



(4.7) 
(4.8) 

(4.9) 



(4.10) 



(4.11) 



(4.12) 
(4.13) 

(4.14) 



equation (4.6) can be represented as one-dimensional integrations as derived in Appendix 
C: 
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1 



87r2i?3 Jo 
e~ 



dxf[x)e ^ ^^evic 



V — 



dxf{x)H,, 



V — ^x 



m — 1 



exp 



[X — ^Vj 

2(1-7^) 



(m = 0) 
(4.15) 

(m > 1) 



which have to be integrated numericaUy. In the high v hmit, these coefficients reduce to 

1 (Ji 



12V37r2 (Jo 



(4.16) 



To see the behavior of two-point correlation function of density peaks on large scales 
which was described in Bardeen et al. (1986), we consider the case, r/j^^) — > 0. The 
two-point correlation function of peaks calculated by our method is, to lowest order in 



(G)/ 
pk \ 



(2), 



(4.17) 



The corresponding expression in Bardeen et al. (1986) [equation (6.12) of their paper] is 

epk=(5)V\ (4.18) 



where 



(5) 



f27r)2i?3n(«)' 



2 / 

dae~" ' I dxf{x) 



exp 



[x — 7a)^ 
2(1-7^) 



a — ^x 



[27r(l - 72)]i/2 1-72- 



(4.19) 



The above two expressions are equivalent. In fact, explicit calculation (see Appendix D) 
shows 



In high f limit, i^pk = v'^tp^'^'^ which corresponds to Kaiser's result (Kaiser 1984). 

Jensen and Szalay (1986) give the expression of the autocorrelation functions of regions 
in which the value of a field a is higher than a threshold v . Autocorrelation functions of 
the thresholded regions are regarded as approximations for correlation functions of peaks 
in high v limit in the context of biased structure formation in astrophysics. Their formula 
can be derived by our method as shown in the previous section: equations (3.19) and 
(3.22) derive their formula [eq. (4) of their paper]. 



Piv(ri, 



Pi 



N 



m = {mi^i\m} k<l 



(4.21) 



r=l 
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where 



/?(sc) 

J. 



(sc) 



we'" 



(m = 0), 
(m > 1). 



(4.22) 



V2erfc(///V2) 

Note that the definition of the Hermite polynomials are different from that in Jensen & 
Szalay (f986). In the case of density peaks, it is obvious that the above quantities A^^^'^ 
should be replaced as follows: 



r I 



^(Pk) 



dxf{x)H^_i i ^ I exp 

\Vl-7/ 



[x — 

'2(1-7^) 



(1 _ ^2^m/2 



dxf[x)e ^^^"^eiic 



V — ^x 



(m = 0), 
(m > I). 

(4.23) 



This quantities are easily evaluated by one-dimensional numerical integrations. This ex- 
pression for A^-point correlation function [eq. (4.21) and eq. (4.23)] of density peaks in 
Gaussian random fields is efficient than the methods developed in Bardeen et al. (1986), 
especially for higher N . In practice, the expansion is summed for arbitrary accuracy as 
in the Jensen and Szalay's formula. Figure 6 shows the various approximations for the 
two-point correlation function of peaks including the approximations of Jensen and Szalay 
[eq. (6) of Jensen & Szalay (1986)], Bardeen et al. [eq. (6.22) of their paper with n = 2], 
peak-background splitting method developed by Bardeen et al. [eq. (6.46) of their paper], 
Lumsden et al. (1989)^, comparing with our method [eq. (4.21) and eq. (4.23)]. In the 
plot, the underlying fiuctuation is assumed to be random Gaussian. We use the smoothed 
power spectrum of standard cold dark matter model given in Bardeen et al. (1986) with 
parameters = I, A = 0, and a Hubble constant of 50 km/s/Mpc. The smoothing length 
R is 0.2/i~^Mpc and the Gaussian window (2.2) is adopted. Figure 7 shows the three- 
point correlation functions of peaks in equilateral configuration for the approximations of 
Jensen and Szalay [eq. (4) of their paper with = 3], peak-background splitting method 
[eq. (6.53) of Bardeen et al. (1986)], Jensen and Szalay with effective thresholds, compar- 
ing with our method for the same underlying spectrum as in Figure 6. In the effective 
threshold method, we use the same effective threshold as in the two-point correlation 
function. In Figure 8 plotted the normalized three-point correlation function 



^ e^'H^l, ^2)e('H^2, rs) + e('H^2, r3)e('H^3, ri) + ^(^H^S, ri)e(2)(ri, 7-2) 

for approximations same as in Figure 7. 



(4.24) 



^They used Jensen and Szalay's formula with effective threshold v^n which is determined by matching 
the two-point correlation function of Jensen and Szalay with that of peaks in large-separation limit, i.e., 
Veff is the solution of the implicit equation, A^l'^\veff) = {a){v) for each true threshold v. 



- 20 - 



4.3 Weighted density extrema of non-Gaussian random field 



The approximation adopted in the previous subsection does not work when the back- 
ground fluctuation is not Gaussian random held. This is because in calculating the biasing 
of the non-Gaussian random held, the higher order correlations at the same point includ- 
ing derivatives such as {a'^ (x) Aa[x)) can not reasonably neglected. In this subsection, 
we consider the density extrema with weighted factor according to the second derivatives 
of the held (see below) above some threshold which can serve as an approximation 
of the density peaks for high threshold. The method to evaluate the density extrema of 
non- Gaussian random held to arbitrary accuracy is given in the following. 

The number density of weighted density extrema above threshold p we consider in this 
subsection is 

p,,t{r) = e{a{r) - p)6^ {d,a{r)) {-if det {d,d,a{r)) , (4.25) 

[cf. eq. (4.4)]. This number density held is for density extrema with weight factor 
— sign[det [didja[r))]. For large threshold the number of maxima dominates the num- 
ber of minima and saddle points, and equation (4.25) is a reasonable approximation for 
density peaks because the weight factor for density maxima is -|-1. Otto et al. (1986) and 
Cline et al. (1987) investigated the average number density and the correlation functions 
of this type of weighted extrema for Gaussian random flelds approximately. We are in the 
place to generalize these analyses for non-Gaussian random flelds. 
The coefflcients for the semi-nonlocal bias (4.25) is 

xe-""!^ [Hk+2{i^)U{p^j]) - Hk+i{v)Ji{{p.,]) 

+ Hk{v)J2{{p^,}) - Hk-l{v)J3{{p^,})] • (4.26) 

The derivation of this result is given in Appendix E. In the above notation, the number 
of indices corresponding to a, da/dx\ d^a/dx^dx^ [i < j) among /ii, . . . , are A;, 4-, p^j, 
respectively. For A; = 0, we use the notation, 

H_,{.)^^e^y^edc(^j=y (4.27) 

Table I gives the deflnition of Jq, Ji, J2 and J3 as functions of {pij}- For cases not listed in 
Table I, J^'s are all vanish. We can calculate the number density and the correlation func- 
tions of weighted density extrema of non- Gaussian held to arbitrary accuracy in principle 
by using the coefflcients (4.26) and rules C-i) ~ C-iv) provided that the diagrammatic ex- 
pansion converges. As primary examples, the number density of weighted extrema nex{i^) 
and the two-point correlation function ^ex,iy{'>^) of Gaussian random flelds to leading order 
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in t/j are, from diagrams in Figure 9, 



2[V 



6 



C(r) 6//2(//2 - 3) AC(r) 



1 



^0^ 



11 



+ 



0-1 



1 



(4.28) 



(4.29) 



where (^(r) = crQt/j[r). Otto et al. (1986) and Catelan et al. (1988) derived the leading 
and sub-leading contribution of p to ^ex,iy Our method provides a way to calculate the 
same quantity not relying on the expansion in p. 

In the high p limit, the leading contribution of p in each diagram is come from the 
term which maximize k for each external point as seen from the expression (4.26). Thus, 
considering the leading contribution of p is equivalent to considering the local bias with 
coefficients 

1 



Rk 



Substituting these coefficients to equation (3.17), we obtain 

N 



P. 



N 



1 



f27r) 



p'e-^' 



(4.30) 



(4.31) 



The number density is 



Pi 



I- 



(71 



2 -v^ 

p e 



Z,{p) 



(4.32) 



which was previously derived by Catelan et al. (1988) by using a different method. The 
correlation function is calculated by taking the connected part of the following quantity. 



P, 



N 



Pi 



N 



Zn{J, = p) 

[Ziwr ' 



(4.33) 



which is the same as in the biasing with sharp clipping in high p limit [eq. (3.25)]. 



4.4 Biasing and hierarchical underlying fluctuations 

In this subsection, we consider the case that the underlying fluctuation satisfy the 
hierarchical model of correlation functions (White 1979; Fry 1984); 

JV-l 

e(^)(a.l,...,a.^.)= ^ gi^) ^ 11 i^'K^A.xs). (4.34) 

trees (a) labelings edges (AB) 
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In the above symbolic notation, the edge (AB) is one of the edges in a tree graph (a) 
which is a set of connected — 1 edges hnking N points, a;i, . . . , without making 
any loop. All the distinct tree graphs are labeled by (a) and "labelings" indicate the 
symmetric sum with respect to the A^-points (Fry 1984). Q^^"^ are constants for each 
tree topology (a). If Q^^"^ = for all A^ > 3, the hierarchical model reduces to random 
Gaussian fluctuation. Fry & Gaztanaga (1993) showed that the local bias which can be 
expanded as Taylor series by density contrast 6 of background field. 



a'. 



Pb{^) = Y.1^^^ (4-35) 

leads the hierarchical form (4.34) of biased field in the leading order of ^^"^^ In this 
argument, the coefficients a'^ are assumed to be the zero-th order of ^^"^^ This argument 
can not be applied to biasing by sharp clipping (3.21) because sharp clipping can not be 
expanded as Taylor series. Moreover, even if we approximate the sharp clipping by some 
smooth function, it is a function of 8 j a and the coefficients a'^ are no longer the zero-th 
order of i^^^^ ~ 0{a^). We complement, in the following. Fry and Gaztanaga's argument 
by investigate the case biasing is a local function of normalized density contrast 8 j a which 
has not to be expanded as a Taylor series as in the sharp clipping bias or in density peaks 
in the approximation of the preceding subsection. 

In the following, the background field a is identified with the normalized density 
contrast 81 a. The field 8 is also assumed to satisfy the hierarchical model (4.34). The 
j-vertex, in this case, contributes the order (t^~^ and we need not to consider any vertex if 
we only consider the leading order in a. Fvaluating diagrams in Figure 10, the following 
results are obtained to the leading order both in a and tp^'^'^: 

Pi = i?o, (4.36) 
&={^)\i2. (4.37) 



eS23 = ^(V'i2V'23 + cyc.), (4.38) 



i?3 

Cii%34 = %?^(V'i2V'23V'34 + sim.(12)) + ( V'i2 V'i3V'i4 + sim.(4)). (4.39) 

The higher orders of r/j^^^ is omitted for we are interested in large-separation limit. Thus, 
the parameters Qjsj I of hierarchical model of a biased field are 



JB) _ R0R2 

.3,« - ^2 ' 

)(B) _ J^O J^2 ^(B) _ ^0^3 

'4,/3 — p4 ' ^4,7 



(4.40) 
(4.41) 



to leading order. The indices of tree graphs a, /3, 7 are defined in Figure 11. For arbitrary 
A^, it is obvious from above calculation that satisfy the hierarchical model to leading 
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order in i^g ~ 0{cr^) and, moreover, the parameters Qjvi t)e directly obtained from 
the topology (a) by simple rules as follows: let nij be the number of external points which 
is attached by j lines in a connected tree graph (a) = ^7 ji^j = 2(A^ — 1)), 

then 

.Ox -2(jv-i) n P™^ 

'2"°" = (I) ^ = n 4". (*■«) 

where = Rj/Rq. In the biasing by sharp clipping, equation (4.22) is applied. In the 
biasing by density peaks, equation (4.23) is applied in the approximation of the preceding 
subsection. In high p limit, both biasing via sharp clipping and via density peaks have 
the same behavior, Ajn{v) ^.^^ p"^ ^ and always Q^^\ = 1. 



4.5 Gravitational evolution as nonlocal biasing 

In our formalism, the gravitational nonlinear evolution of the density field can be seen 
as a kind of nonlocal biasing from the initial fiuctuation field. In this subsection, we 
give explicit relation between our formalism and the perturbation theory of gravitational 
instability and apply our formalism to the calculation of correlation functions of density 
field. We consider the non-relativistic coUisionless self-gravitating system in the fiuid 
limit for Einstein-de Sitter universe (0 = 1, A = 0, p = 0) as an example. The mass 
density contrast 8{x,t) and the peculiar velocity field v{x,t) are governed by the equations 
(Peebles 1980): 

^ + lv.[(l + ^)^] = 0, (4.43) 
01 a 

dv a 1 , 1 „ , , 

— + -v + -{vV)v + -V(i> = {), 4.44 
Ot a a a 

A(f) = 4:TrGpaH, (4.45) 

where V denotes the derivative by means of co moving coordinate a;, p and a[t) are the 
average density and expansion factor, respectively, and are given hy p = (GTrCt^)"^, 
a oc t^/"^ in the Einstein-de Sitter universe. The linear solution is derived by linearizing 
equations (4.43)-(4.45). The growing mode of the linear solution has the form, 

6('\x,t)=^^Y\{x), (4.46) 

where tj. is initial time or recombination time, e(a;) is interpreted as the primordial fiuc- 
tuation. The second order perturbative solution for growing mode (Peebles 1980) is given 
by 

6('\x,t) = (^^y [^e'{x) + Ve{x) • VU{x) + ^ [d,d,U{x)] mU{x)]^ , (4.47) 
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where 



The mass density, 



Uix] 



1 

47r 



\x — X' 



p^,{x,t) = -p{t)[l + 6^'\x,t) + 6^'\x,t) + 0{e% 



(4.48) 



(4.49) 



is regarded as nonlocal biasing where the background field a is identified with normalized 
primordial fiuctuation e/cr, where cr^ = (e^). The kernels are derived to be 



A^(°) = Pit), 

/ / \ 2/3 

K(^\x) = a i^-^j pit)S^ix) + Oia^) 



/ f \ 4/3 

K^'H^,y) = <^\^) Pit) 



A-W = 0(a") (n > 3), 



(4.50) 
(4.51) 

1 dS''{x)d\y\-^ 1 9(53(2/) a 1*1 -1 



iir dyi dxi 



+ 



1 d^\x\-^ d^\y\-^ 
287r2 dxidxj dyidyj 



(4.52) 
(4.53) 



which are nonlocal. Using these kernels, and assuming that the primordial fiuctuation e 
is Gaussian, our rules A-i) ~ A-vi) for two- and three-point correlation function give the 
following results after lengthy algebra: 



Lt{ri,r2;t) 



Cpt{ri,r2,r3;t) 



-\- cos 02 



6n(ri2) + Oia') 



(4.54) 



7 



(5 + 2cOs2 ^2)6n(ri2)6n(r23) 



Uri2)Jt{r2s) CL{r2s)Jnri2] 
9 ~r 9 



23 



12 



,4/ .. 4"(^i2)Jjr(r23) 
+ - 3(3 cos ^2 - 1) -^-^ 

I ' 1 9' 9 



12' 23 



+ (1 -3cos2 



6n(ri2)4"(^23) 6n(r23)Jjr(^12) 



23 



12 



+cyc.(3 terms) + 0{a^). 
In the above, the two-point correlation function of the primordial field is 

6n(ri2) = (e(ri)e(r2)) 



(4.55) 



(4.56) 



where ri2 = l^'i — T2I; 62 is an angle between two vectors, Vi — 7*2 and 7*3 — 7*2; and 

Tin 
'3 



J™(r) = Jo dr'^^ini'i''). Recently, Bharadwaj (1994) independently obtained the same 
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result for three-point correlation function (4.55). Putting oc r"''" in equation (4.55), 
we rederive Fry's result (Fry 1984). For equilateral configuration, ri2 = = rsi, the 
hierarchical amplitude, Qpt = Cpt/(3(^pt) is 

Qpt (equilateral) = ^^'^ ■ (4.57) 

'{'i - 7inj 



5 CONCLUSIONS 



We have developed a formalism to explore the statistics of biased field in large-scale 
structure of the universe. In particular, the general nonlocal biasing can be treated in this 
formalism. The underlying field is not necessarily assumed to be a random Gaussian field. 
The formalism is based on diagrammatic series expansion of the correlation functions of a 
biased field. As for local biasing, the form of series expansion of a biased field is known by 
previous works, but it is very complicated expression (Borgani & Bonometto 1989, 1990). 
Our formalism, when applied to local biasing problem, simplify this complexity. 

In quantum field theory and statistical physics, the Feynman diagrams play roles 
more than mere simplification of calculations. It would be interesting if our formalism 
could play similar roles. In fact, our diagrammatic expansion can be conceived to be the 
generalization of the Fdgeworth expansion which is useful in connecting dynamics and 
statistics in large-scale structure by the argument in §4.1. 

We have applied our formalism to the problems which have been studied previously 
to show the effectiveness of our method. Application to peak statistics improves the tech- 
nique developed by Bardeen et al. (1986). In particular, we can evaluate the higher order 
correlations of peaks successfully for which technique of Bardeen et al. had difficulty. For 
statistics of weighted extrema, which could be an approximation to peak statistics, we 
showed the general way of calculation not depending on the high-threshold expansion. 
Weighted extrema is observable in surveys of large-scale structure and our result can be 
directly compared with observations. The similar technique was recently applied to the 
genus statistics of weakly non-Gaussian field by the present author (Matsubara 1994). 
The conservation of hierarchical model of correlation functions subjected to local biasing, 
which is studied by Fry and Gaztanaga (1993), was revisited using our formalism. We 
complement their results by investigating a local biasing through normalized density con- 
trast which may not be represented by Taylor series, as in sharp- clipping biasing. The 
gravitational evolution of primordial fiuctuation also falls under the category of nonlocal 
biasing as shown in §4.5. 

Those problems have been studied individually before this paper. Our formalism 
enables more detailed examination of those problems. Those examples show that our 
formalism is not only a general method but also a practically powerful way to investigate 
various problems. The formalism developed in this paper may be an infiuential method 
for future investigations of large-scale structure of the universe. 



- 26 - 



ACKNOWLEDGEMENTS 



I thank Y. Suto, and J. Yokoyama for helpful discussions. I thank T. Watanabe for 
his help in plotting Figure 6, 7 and 8 and discussions. 

I acknowledge the support of JSPS Fellowship. This research was supported in part by 
the Grants-in-Aid for Scientific research from Ministry of Fducation, Science and Culture 
of Japan (No. 0042). 



- 27 - 



APPENDIX A: EXPECTATION VALUES OF PRODUCTS OF 



GENERALIZED WIENER-HERMITE FUNCTIONALS 

Imamura, Meecham & Siegel (1965) gives the method to calculate symbolically the 
expectation values of products of Wiener-Hermite functionals. In this Appendix, we gen- 
eralize their method to the case of generalized Wiener-Hermite functionals. The quantity 
we need to calculate is 

\t=l I G 

To derive the method to calculate this quantity, there are at least two ways, i.e., using 
operators and using generating functionals. We present both ways here although either 
is sufficient to our purpose. 

Al The method using operators 

The method using operators is to use the operator representation of generalized 
Wiener-Hermite functionals. This method is similar to the derivation of Feynman rules 
in operator formalism in quantum field theory. 

We define the annihilation operator a{x) and creation operator a'^{x) as 

= \ I d^y^l^~\x, y)a{y) + (A.2) 
z J 6a[x) 

a\x) = la{x) - [ d^yij{x,y)-^ (A.3) 

Annihilation and creation operators transform as covariant and contravariant vectors, 
respectively under the linear transformation (2.19). These operators indeed satisfy the 
following commutation relations as annihilation and creation operators: 

[a{x),a\y)] =8\x-y), (A.4) 
Hx),a{y)] = \a\x),a\y)\ =0. (A.5) 

In terms of annihilation and creation operators, the generalized Wiener-Hermite function- 
als are represented by 

m 

7i(„) (a?!,. ..,«„) = : J]^ [a(a;j) a^(a;j)] :, (A. 6) 

8 = 1 

where 

a{x) = j d^yil){x,y)a{y), (A.7) 
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is a contravariant vector. The notation : • • • : means the Wick product (Wick 1950): place 
all the annihilation operators to the left of all the creation operators in the Wick product. 
For example, 

: \a{x) + a\x)\ \a{y) + a.\y)\ '. = a{x)a{y) + a{x)a\y) + a{y)a\x) + a\x)a\y), (A. 8) 



and so on. The representation (A. 6) is proved by induction using the recursion relation 

(2.23) . 

We then define the "vacuum state" |0) as 

\0) ^ {PG[a]V^\ (A.9) 

where Paict] is an infinite dimensional Gaussian distribution functional given by equation 

(2.24) . The "n particle states" (n > 0) are defined using vacuum state as 



1 ^ * * • •) *^ ' 



= a^{x^)---a\x^)\{)), (A.IO) 



which are the bases of the Fock space JF. The following properties hold for n particle 
state: 

a^y)\xi,...,Xn) = (A. 11) 

a(j/)|0)=0, (A.12) 

n 

a{y)\xi, ...,Xn) = J2^^iy ~ • • • • • (A. 13) 

8 = 1 

We define the inner product of n particle state \xi^ . . . , Xn) and m particle state . . . , y^) 
in the Fock space by 



(a?!, a;„|?/i, ...,?/„) = 4„ J d^z^ ■ ■ ■ d^z.^^) ^{x^, z^) ■ ■ ■ 4) ^(a;„, z„) 

X y"[Ja]|zi,...,z„)|j/i,...,j/„), (A. 14) 

where we employ Dirac's notation for the inner product. The annihilation operator and 
the creation operator are Hermite conjugate to each other under this definition of the 
inner product in the Fock space. The inner product (A. 14) can be calculated explicitly 
by commutation relations (A. 4), (A. 5) and the properties (A. 13) resulting in 



(«!,..., a;„|j/i,...,j/„) = (0 



a{xi) ■ ■ ■a{x.a)a\y-^) ■ ■■a\y^ 







= [(53(a;i - J/i) • • • <5^(a;„ - y^) + sym.(j/i, . . . , y^)\ , (A.15) 

where sym.(j/;^, . . . , y,^) represents the symmetric summation of the previous term with 
respect to y^^ . . . , y^. If F[a,a^] is a c-number, i.e., the operator which simply multiplies 
the real number to the operand, then 

(F[a,at])^ = (o|F[a,at]|o). (A.16) 
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The quantity (4.53) is, after all, represented as 



\j=i 

mi 



8 = 1 8 = 1 



0) . (A.17) 



Expanding the products in r.h.s, and using the commutation relations (A. 4), (A. 5) and the 
properties (A.11)-(A.13), this quantity is calculated algebraically. But more convenient 
diagrammatic calculation can be introduced. The "particle" created by a creation operator 
should be annihilated by an annihilation operator placed left of that creation operator. 
In a single Ti-(m,) represented by Wick product, the creation operators are already placed 
to the left of annihilation operators, and the "particles" created in some Ti-(m,) can not 
annihilated by the annihilation operator in the same ?i(m,)- These observations show the 
following diagrammatic rules for calculating equation (A.17): 

i) Corresponding each draw rrij points labelled by x^j^\ . . . , x^j^^\ 

ii) Make ^jf^ pairs out of those points such that the two points in the same Ti-(mj) 
are not paired. If ^jf^ odd number, the equation (A.17) vanishes. 

iii) Associate factor i/j^x^J'^ ^i'^) foi" each pair, x^J''^ and a;|,'^ and make products of these 
factors. 

iv) Sum up those products obtained from all the possible pairings. 
For example, the quantities. 



^(i)(a^)^(i)(?/)^(2)(^,«^))^, (A.18) 

are evaluated by the diagrams in Figure 12. Applying the above rules we obtain 

V'(a;,z)V'(j/,«?) + tl^{x,wyi;{y,z). (A. 19) 

When the generalized Wiener-Hermite functionals of type ?-{(™) are contained in the l.h.s. 
of equation (A.17), the factor iIj{x''J'\ «i'^) or (5"^(a;*j'') — «i'^) or i/j'^^x'^J'^ «i'^) are associated 
for each pairs and which factor should be associated can be determined by the fact that 
the result should be transformed correctly under linear transformation (2.19). 
The following orthogonality relation can be derived by the above rules: 



n^n){xi, . . . , Xn)n^''^\y^, . . . , 



G 



8.mn [^'^{xi - J/i) • • • <5^(a;„ - y^) + sym.(j/i, . . . , y^)] . (A.20) 
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A2 The method using generating functional 

The second method using generating functionals is to derive rules i) ~ iv) using gen- 
erating functional of generalized Wiener-Hermite functionals. This method is similar to 
the derivation of Feynman rules in path integral formalism in quantum field theory. 

The generating functional of generalized Wiener-Hermite functionals is defined by 



so 



that 



OO 1 . 

rTn rnl J 



m = 



\Xi^ . . . , X,f 



J=0 



8J{xi) ■ ■ ■ 6J[x 

Lctional: 

F[J -a]= exp {-^ / d^xd^y[J{x) - a{x)]ij-^{x, y)[J{y) - a{y)] 



We consider the following functional: 

1 



and expand this functional with respect to J: 

F[J - a] = £ -L / d^., . . . d^^Jix,) . . . Jix^) .j^^'Tsj^.^ 



m = 



ml 



(A.21: 



(A.22) 



(A.23) 



(A.24) 



j=o 



Because 



a 



6J{xi) ■ ■ ■ 6J{x, 



j=o 



6a{x-i 



a 



6a{xi) ■ ■ ■ 6a{xj 



■ Saixr 



F[a\H^ 



j=o 



m). 



X 



(A.25) 



j=o 



we see 



F[J 



F[a] 



= exp J d^xd^yJ{x)ij-^x,y)J{y) + J d^xd^ya{x)ij-^x,y)J{y)^ (A.26) 
Using this representation of generating functional, we can derive the following equation: 

J2 J d^xd^yJ,{xyi;-\x,y)J,{y) 



{Q[a, Ji] • • • Q[a, Js])q = exp 



The rules i) ~ iv) are derived by operating 



(A.27) 



8,h {x^^^) ■■■8,h {x^r^) SJ, ...SJ, {xi-'^) ' 

to the both side of equation (A.27), and putting Ji = J2 = • • • = = 0. 



(A.28) 
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APPENDIX B: THE CALCULATION OF THE STATISTICAL FACTOR 



In the expression (3.5) the points asso dated with 'H^''^+-+''^\'H(^m, ) are classified into 
m -\- N groups with ni, . . . , n^; ?7^i, • • • , rn]\j members respectively. These points in each 
group are symmetric about permutations. We distinguish these m-\-N groups by labels like 
a, /3, and so on. In evaluating the Gaussian average of the product of generalized Wiener- 
Hermite functionals using diagrammatic method described in Appendix A, we denote the 
number of lines connecting the members of two groups a and /3 as a symmetric matrix 
Maf5 with vanishing diagonal elements. Obviously, graphs which have the same value of 
the matrix Af„^ give the same contribution to Pjv- Thus we need to calculate the degree 
of multiplicity of the graphs which have the same value of Af„^. 

We designate the number of points in a group a as Af„. The number of ways to divide 
the Af„ points into Af„^i, . . . , Af„^m_|_jv parts for all groups a's is 




The number of ways to connect these points is given by na</3 ^a/?! and thus the multi- 
plicity of the graphs is 

• • • n„!mi! • • • mjv! 
ff • (^-2) 

J_ J_ ^ends 
ends 

The numerator cancels out in equation (3.5) and the statistical factor (3.6) follows. It is 
obvious that the m vertices are distinguished in the above calculations and the graphs 
like Figure 13 are counted as different graphs. 



APPENDIX C: EVALUATION OF 



This Appendix is devoted to the derivation of equation (4.15). The following identity 
is convenient for our purpose: 



2 2(1 -72) ~ 2 2(1 -72)' 



(c.r 



From this identity, equation (4.8) reduces to 



1 



f27r)2 R3 7o 



dxf[x)e' 





[p — 7a;)2 


exp 


[ 2(1-72)J 



[2vr(l - 7^)]^/^ 



(C.2) 
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and the integral in equation (4.7) can easily be calculated: 



1 



87r2 i?3 Jo 



dxf[x)e ^ ''^erfc 



V — 

2(1-7^), 



(C.3) 



The latter expression (C.3) was previously derived by Mann, Heavens & Peecock (1993). 
Differentiating equation (C.2) and using the identity (C.l) again, we obtain 



dv 



m — l 



(27r)5/2i?3(l _ ^2)m/2 70 

for m > 1. Thus, 



dxf{x)H^_i ( ^ ) exp 



[X — ^Vj 

2(1-7^) 



,(C.4) 



1 

87r2i?3 7o 

e 



-1/2/2 



I (27r)5/2i?3(l _ ^2)m/2 70 

is proved. 



V2(l-7^)/ 
dxf{x)Hm-i 1^-^^==^^ exp 



2(1-7^) 



(m = 0) 
(C.5) 

(m > 1) 



APPENDIX D: PROOF OF 



(d>=AA,:«V)/n!i:V 



pk 



In this Appendix, we calculate (a) defined by equation (4.19). Using the identity 
(C.l), the equation (4.19) reduces to 



(5) 



2 / 

dxf{x)e~'' / da 



exp 



[a — 7^)^ 



2(1-7^) 



a — ^x 



[27r(l - 72)]i/2 



fD.r 



The integration with respect to a is straightforward by a transformation jS = a — ^x. 
Comparing the result with equation (C.2), this quantity has the simple form. 



(«) = (G), . 



fD.2l 



npk I'^J 



- 33 - 



APPENDIX E: DERIVATION OF COEFFICIENTS FOR DENSITY 

EXTREMA 



In this Appendix, we derive the coefficients (4.26). The quantity we should calculate 



IS 



Rik,{k},{p,,}) 



n^aftj n^a|-)'"''(«-)«'(ft)(-i)'det,K,)) (E.1) 



1 



(27r)VdetM 



da ]^ dfii ]^ diOij exp 



8 = 1 i<j 



1 



where ten dimensional vector A,, is 



'''*(°-'''n^nf^V"det.. 



(A^) = (a,/3i,/32,/33,^ll,^22,^33,^23,^13,^12)- 

The correlation matrix M^^ = {A^A^) is explicitly 



M 



1 


O^v 




ol 1 






o 


o 




o 


15(7„2 


o 


Ov 


o 


o 


15^2^ J 



where 





" " 




" 1 " 




" 





" 




" 1 





" 




CO 


1 


1 " 


Oy = 





= 


1 


,o = 
















1 







1 


3 


1 









1 





















1 




_ 1 


1 


3 



(E.2) 



(E.3) 



(E.4) 



(E.5) 



This is the consequence of spatial homogeneity and isotropy (see Bardeen et al. 1986). 
We introduce new variables cu,-,- as 



(J, 



then the non- vanishing correlations among variables a, are 



(E.6) 



(7^ 



2 ' 
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= = = (E.7) 

and all other correlations vanish. The coefficients (E.l) is represented by new variables, 

a, (3,, CO, J as 



where 



+ --^ (^23 + ^1^3 + ^1^2 - ^22^33 " ^11^^33 " ^11^22) « 
+tJlltJ22^33 + 2tJ23^13^12 " ^^11^2^3 " ^22^^13 " ^33^lV (E-9) 

We see, from equation (E.9), that the equation (E.8) survives only when integers pij take 
the values listed in Table 1. Using Ji{{pij}) defined by Table 1, the multivariate Gaussian 
integral by uj of the second factor of the equation (E.8) reduces to 



.2 \ 3-J2,<jP'3 I f)k 



a — V \ 



a 



[H^{a)M{p,,]) - i?2(«)Ji(te,}) + 



+H,{a),M{p,,}) - ,M{p.,})] . (E.IO) 



/ G 



Defining 



i7_,H = y|e''VVfc(-^), (E.ll) 



we 



find 



'd''9{a-u) „ , A 1 



\ dak -V ^ 
for n = 0, 1,2, . . .. The first factor of equation (E.8) is 



H^{a)) = —Hk+n-iiu)e-^ ^\ (E.12) 



^ ' • HiMHiMHisiO). (E.13) 



(27r)3/2 VV3^ 
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where 

^'<°)={?-i,'/^(,-i)!! I'™!;. ■ 

Therefore, we finally get the explicit form of the kernels: 

xe-""!^ [Hk+2{i^)U{m]) - Hk+i{v)Ji{{p.,]) 

+Hk{v)J2{{p^,}) - Hk-^{v)U{p.,})] . (E.15) 
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Pll 
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P33 


P23 


Pl3 
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n 
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9 


n 


n 


n 


n 


9 


n 
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-2 




















2 








-2 





1 


1 


1 




















1 


1 








2 

















-2 





1 








2 














-2 








1 








2 











-2 











1 


1 


1 











2 



Table 1: The definition of .]{. For other values of not listed in this table, .]{ = 0. 



- 37 - 



REFERENCES 

Adler R. J., 1981, The Geometry of Random Fields (Wiley: Chichester) 
AUen T. J., Grinstein B., Wise M. B., 1987, Phys. Lett., B197, 66 

Appel P., de Feriet J. K., 1926, Fonctions Hypergeometriques et Hyperspheriques, Polynomes 

d'Hermite^ Ganthier-Villars, Paris 
Babul A., White S. D. M., 1991, MNRAS, 253, 31P 
Bahcall N. A., Soneira R. M., 1983, ApJ, 270, 20 

Bardeen J. M., Bond J. R., Kaiser N., Szalay A. S., 1986, ApJ, 304, 15 

Bardeen J. M., Steinhardt P. J., Turner M. S., 1983, Phys. Rev. D, 28, 679 

Bernardeau F., Kofman L., 1994, preprint IFA-94-19 (astro-ph/9403028). 

Bertschinger F., 1992, in New Insights into the Universe^ eds. Martinez, V. J., Portilla, M. 

& Saez, D., Springer- Verlag, Berlin-Heidelberg 
Bharadwaj S., 1994, ApJ, 428, 419 
Borgani S., Bonometto S. A., 1989, A&A, 215, 17 
Borgani S., Bonometto S. A., 1990, ApJ, 348, 398 

Bower R. G., Coles P., Frenk C. S., White S. D. M., 1993, ApJ, 405, 403 
Catelan P., Lucchin F., Matarrese S., 1988, Phys. Rev. Lett., 61, 267 
Cen R., Ostriker J. P., 1992, ApJ, 399, L113 
Cen R., Ostriker J. P., 1993, ApJ, 417, 415 

Cline J. M., Politzer H. D., Rey S.-J., Wise M. B., 1987, Commun. Math. Phys., 112, 217 
Coles P., 1989, MNRAS, 238, 319 

Davis M., Ffstathiou G., Frenk C. S., White S. D. M., 1985, ApJ, 292, 371 
Davis M., Peebles P. J. F., 1983, ApJ, 267, 465 

Doroshkevich A. G., 1970, Astrophysics, 6, 320 (transl. from Astrofizika, 6, 581 [1970]) 
Ffstathiou G., Kaiser N., Saunders W., Lawrence A., Rowan- Robinson M., Fills R. S., Frenk 

C. S., 1990, MNRAS, 247, lOP 
Fry J. N., 1984, ApJ, 279, 499 
Fry J. N., 1986, ApJ, 308, L71 
Fry J. N., Gaztanaga F., 1993, ApJ, 413, 447 
Grinstein B., Wise M. B., 1986, ApJ, 310, 19 
Guth A., Pi S.-Y., 1982, Phys. Rev. Lett., 49, 1110 
Hauser M. G., Peebles P. J. F., 1973, ApJ, 185, 757 
Hawking S., 1982, Phys. Lett., 115B, 295 

Hodges H. M., Blumenthal G. R., Kofman L. A., Primack J. R., 1990, Nucl. Phys., B335, 
197 

Ikeuchi S., 1981, PASJ, 33, 211 

Imamura T., Meecham W. C, Siegel A., 1965, J. Math. Phys., 6, 695 



- 38 - 



Jensen L. G., Szalay A. S., 1986, ApJ, 305, L5 

Juszkiewicz R., Weinberg D. H., Amsterdamski P., Chodorowski M., Bouchet F. R., 1994, 

ApJ in press (astro-ph/9308012). 
Kaiser N., 1984, ApJ, 284, L9 

Klypin A. A., Kopylov A. I., 1983, Sov. Astron. Lett., 9, 41 
Kofman L. A., Linde A. D., 1985, NucL Phys., B282, 555 

Kolb E. M., Turner M. S., 1990, The Early Universe (Addison- Wesley: New York) 

Lumsden S. L., Heavens A. F., Peacock J. A., 1989, MNRAS, 238, 293 

Ma S.-K., 1985, Statistical Mechanics^ World Scientific, Philadelphia 

Maddox S. J., Ffstathiou G., Sutherland W. J., Loveday J., 1990, MNRAS, 242, 43P 

Mann R. G., Heavens A. F., Peacock J. A., 1993, MNRAS, 263, 798 

Matarrese S., Lucchin F., Bonometto S. A., 1986, ApJ, 310, L21 

Matsubara T., 1994, ApJ, 434, L43 

Moore B., Frenk C. S., Weinberg D. H., Saunders W., Lawrence A., FUis R. S., Kaiser N., 

Ffstathiou G., Rowan-Robinson M., 1992, MNRAS, 256, 477 
Ortolan A., Lucchin F., Matarrese S., 1989, Phys. Rev. D, 40, 290 
Ostriker J. P., Cowie L. L., 1981, ApJ, 284, L9 

Otto S., Politzer H. D., Wise M. B., 1986, Phys. Rev. Lett., 56, 1878 

Peebles P. J. F., 1980, The Large-Scale Structure of the Universe^ Princeton University Press, 
Princeton 

Politzer H. D., Wise M. B., 1984, ApJ, 285, LI 

Postman M., Huchra J. P., Geller M. J., 1992, ApJ, 384, 404 

Salopek D. S., 1992, Phys. Rev. D, 45, 1139 

Salopek D. S., Bond J. R., 1990, Phys. Rev. D, 42, 3936 

Salopek D. S., Bond J. R., Ffstathiou G., 1989, Phys. Rev. D, 40, 1753 

Saunders W., Rowan- Robinson M., Lawrence A., 1992, MNRAS, 258, 134 

Scherrer R. J., Bertschinger F., 1991, ApJ, 381, 349 

Starobinskh A. A., 1982, Phys. Lett., 117B, 175 

Szalay A. S., 1988, ApJ, 333, 21 

Vilenkin A., 1985, Phys. Rep., 121, 263 

White S. D. M., 1979, MNRAS, 186, 145 

Wick G. C., 1950, Phys. Rev., 80, 268 



- 39 - 



FIGURE CAPTIONS 



Figure 1 : Diagrammatic rules for the nonlocal bias. The open circles and filled circles 
corresponds to external points and vertices, respectively. 

Figure 2 : Diagrammatic rules for the local bias. 

Figure 3 : Diagrammatic rules for the semi-local bias. 

Figure 4 : Diagrammatic rules in Fourier space. 

Figure 5 : Diagrams for the Fdgeworth expansion. 

Figure 6 : Various approximations for the two-point correlation function of peaks in CDM 
model [h = 0.5) for thresholds p = 1.5^ p = 2.0, p = 2.5, and p = 3.0. The smoothing 
length is 0.2/i~^Mpc and the Gaussian window function is adopted. Solid lines: our 
method. Dotted lines: Jensen and Szalay (1986). Short Dash lines: peak-background 
splitting. Long dash lines: Bardeen et al. (1986). Dot-short dash lines: Lumsden 
et al. (1989). Dot-long dash lines: original correlation function of CDM model for 
reference. Normalization is not relevant for our analysis. 

Figure 7 : Various approximations for the three-point correlation function of peaks in CDM 
model. The configurations of three-points is equilateral and the length of edges of the 
equilateral triangle is the horizontal axis. Solid lines: our method. Dotted lines: Jensen 
and Szalay (1986). Short Dash lines: peak-background splitting. Dot-short dash lines: 
Jensen and Szalay with effective thresholds. 

Figure 8 : The normalized three-point correlation function of equilateral configuration for 
the same approximations as in Figure 7. 

Figure 9 : Diagrams for n^x. and ^ex,iy 

Figure 10 : Diagrams for hierarchical underlying fiuctuation. 
Figure 11 : Definitions for tree graphs a, /3, 7. 

Figure 12 : An example of diagrams for expectation values of products of generalized Wiener- 
Hermite functionals. 

Figure 13 : These two graphs are counted separately for the statistical factor described in 
the text. 
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